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Abstract. This paper presents the results of a detailed X-ray spectral analysis of a sample of 123 X-ray sources detected with 
XMM-Newton in the Lockman Hole field. This is the deepest observation carried out with XMM-Newton with more that 600 
ksec of good EPIC-pn data. We have spectra with good signal to noise (>500 source counts) for all objects down to 0.2-12 
^ ■ keV fluxes of ~ 5 x lO"'"" ergcm"^ s"' (flux limit of ~ 6 x 10"'^ ergcm"^ s"' in the 0.5-2 and 2-10 keV bands). At the time 

OO ' of the analysis, we had optical spectroscopic identifications for 60% of the sources, 46 being optical type-1 AGN and 28 

optical type-2 AGN. Using a single power law model our sources' average spectral slope hardens at faint 0.5-2 keV fluxes 
' but not at faint 2-10 keV fluxes. We have been able to explain this effect in terms of an increase in X-ray absorption at faint 

fluxes. We did not find in our data any evidence for the existence of a population of faint intrinsically harder sources. The 
average spectral slope of our sources is ~1.9, with an intrinsic dispersion of ~0.28. We detected X-ray absorption (F-test 
significance >95%) in 37% of the sources, ~I0% in type-1 AGN (rest-frame Nh ~ 1.6 x 10^' - 1.2 x 10^^ cm"-) and ~ 77% 
(rest-frame Nh ~ 1.5 x 10^' - 4 x 10^^^ cm"^) in type-2 AGN. Using X-ray fluxes corrected for absorption, the fraction of 
• absorbed objects and the absorbing column density distribution did not vary with X-ray flux. Our type- 1 and type-2 AGN 

do not appear to have different continuum shapes, but the distribution of intrinsic (rest-frame) absorbing column densities is 
different among both classes. A significant fraction of our type-2 AGN (5 out of 28) were found to display no substantial 
5_j ' absorption (Nh < 10^' cm"^). We discuss possible interpretations to this in terms of Compton-thick AGN and intrinsic Broad 

I Line Region properties. An emission line compatible with Fe Kor was detected in 8 sources (I type-1 AGN, 5 type-2 AGN and 

. 2 unidentified) with rest frame equivalent widths 120-1000 eV. However weak broad components can be easily missed in other 

sources by the relatively noisy data. The AGN continuum or intrinsic absorption did not depend on X-ray luminosity and/or 
redshift. Soft excess emission was detected in 18 objects, but only in 9 (including 4 type-1 AGN and 4 type-2 AGN) could we 
■ fit this spectral component with a black body model. The measured 0.5-2 keV luminosities of the fitted black body were not 

' significantly different in type- 1 and type-2 AGN, although the temperatures of the black body were slightly higher in type-2 

AGN «kT>=0.26 + 0.08) than in type-1 AGN «kT>=0.09 ± O.OI). For 9 sources (including I type-I AGN and 3 type-2 AGN) 
a scattering model provided a better fit of the soft excess emission. We found that the integrated contribution from our sources 
to the X-ray background in the 2-7 keV band is softer (F = 1.5 - 1.6) than the background itself, implying that fainter sources 
need to be more absorbed. 

Key words. X-rays: general. X-rays; difi^use background, surveys, galaxies: active 



1. Introduction 

The extragalactic X-ray background (XRB) at energies above 
~ 0.2 keV is made up of the integrated emission from point 
sources, mostly Active Galactic Nuclei (AGN). Synthesis mod- 
els of the XRB (e.g. Setti & Woltjer, 1989; Gilli et al., 20011 
Ueda et al., 12003 ). based on unification schemes of AGN, can 

Send offprint requests to: S. Mateos, e-mail: sm279@star .le.ac.uk 
* Based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA 



reproduce the spectral shape of the XRB with the superposi- 
tion of a mixture of absorbed and unabsorbed AGN. There is 
much observational evidence supporting the unified model of 
AGN (Antonucci, 1993). For example the discovery of large 
columns of X-ray absorbing gas in type-2 AGN (Awaki et al., 
119911 Risaliti et al., 1999j, and the lack of this absorbing ma- 
terial in type-1 AGN. 

With the launch of the Chandra and XMM-Newton obser- 
vatories, our knowledge of the nature and cosmic evolution of 
AGN has increased significantly. However the amount of ob- 
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servational results that cannot be explained in terms of the uni- 
fied model of AGN is also significant. There is a substantial 
number of type-1 AGN for which X-ray absorption has been 
detected (Mittaz et al.. il999. Fiore et al.,^01 Page et al., 
UUUTl Schartel et aL EUHTl Tozzi et al. flUUTl Mainieri et al., 
PWS; Brusa et al.,'2(503 , Page et al.,'M5T, Carrera et al., 2004, 
Perola et al., 2004; Mateos et al.. 2005L as well as Seyfert 2 
galaxies unabsorbed in X-rays (Pappa et al.. l200l| Panessa et 
al.,,2002 Barcons et al.. l20U3l Mateos et al., 2005). The ori- 
gin of the X-ray absorption in type-1 AGN is still not clear. 
Possible explanations include cold gas in the host galaxy. 

For all these objects their optical and X-ray properties can- 
not be explained in terms of an orientation effect only. 

In order to gain insight into these problems and to un- 
derstand better the X-ray emission and cosmic evolution of 
the AGN that make up most of the XRB, spectral analysis of 
large samples of AGN detected in medium and deep X-ray sur- 
veys have been or are being conducted (Mainieri et al., 2002J| 
Piconcelli et al., l2002l l2003l Georgantopoulos et al., .20041 
Caccianiga et al., 20041 Perola et al., 120041 Delia Ceca et al., 
12004 Mateos et al., 20051 However, at the moment, only a 
small number of these studies have carried out a proper spectral 
analysis of the spectra of each individual source, and in many 
cases some assumptions had to be made prior to the spectral 
analysis (frequently on the spectral slope). These studies pro- 
vide observational constrains with large uncertainties. 

The Lockman Hole field is one of the sky regions best stud- 
ied at X-ray wavelengths, because the Galactic absorbing col- 
umn density in this direction is minimal (5.7 x lO'^cm"^, see 
Lockman et al.. ll986t . 

XMM-Newton has carried out its deepest observation in 
the direction of the Lockman Hole field. These observa- 
tions have allowed us to extract good quality (>500 0.2- 
12 keV counts) X-ray spectra for objects down to 0.2-12 
keV fluxes of ~ 5 x 10"'^ ergcm^^sec"' (the flux limit is 
~ 6 X 10"'^ erg cm"^ sec"' in the 0.5-2 and 2-10 keV energy 
bands). We have used the XMM-Newton observations in the 
Lockman Hole to carry out a detailed analysis of the X-ray 
emission of the 123 brightest objects detected in the field. The 
results from the analysis of a sample of fainter objects will be 
described in a forthcoming paper 

Using a subset of these observations, Hasinger et al. (200 IJl 
presented the source detection and properties of X-ray sources. 
Mainieri et al. J2002t conducted a X-ray spectral analysis of the 
objects detected in the field. Worsley et al. (2004 '2005') used 
the total observation of the field to calculate the fraction of un- 
resolved XRB in different energy bands. Finally Streblyanska 
et al. J200 5) have conducted a detailed study of the Fe Ka emis- 
sion in the stacked spectra of type-1 and type-2 AGN. They 
found in their analysis indications for broad relativistic lines in 
both type-1 and type-2 AGN. 

This paper is organised as follows: Sec.|2ldescribes the X- 
ray data that we used for the analysis; Sec. |3l describes how 
we built our sample of objects; Sec. |3 explains how the time 
averaged spectra were extracted for each individual object; in 
Sec.|5lwe show the current status of the optical identification 
process; Sec.|6ldescribes the models that we used to fit the X- 
ray emission of our sources; we show the results of the analysis 



in Sec.0 Sec.|S]and Sec.|5J the dependence of spectral param- 
eters with luminosity and redshift is show in Sec. 1101 In Sec. II II 
we discuss possible explanations for the lack of X-ray absorp- 
tion signatures found in the spectra of five of the sources in our 
sample of type-2 AGN; In Sec. El we compare the integrated 
emission of our sources with the cosmic X-ray background in 
the 2-7 keV band; the results of our analysis are summarised in 
Sec.ini 

Throughout this paper we have adopted the WMAP derived 
cosmology with Hq = 70kms"' Mpc"', Q.m = 0.3 and Qa - 
0.7. All errors are computed with a delta chi-square of 2.706, 
equivalent to 90% confidence region for a single parameter. 

2. XMM-Newton observations 

The deepest XMM-Newton observation has been carried out in 
the direction of the Lockman Hole field, that is centred on the 
sky position RA: 10:52:43 and DEC: +57:28:48 (J2000). The 
data were obtained by adding 17 XMM-Newton observations 
obtained from 2000 to 2002. A summary of the observations 
used in this analysis is given in Table[l] The Lockman Hole was 
also observed during revolution 07 1 , however at the time of this 
analysis there was no Observation Data File (ODF) available, 
and hence, we could not reprocess the data. The first column 
in Tabled shows the revolution number and observation iden- 
tifier. The second column shows the phase of the observation 
(i.e. PV for observations during the EPIC-Payload Verification 
Phase, and AOl and A02 for observations during the first and 
second Announcement of Opportunity). The third and fourth 
columns list the coordinates of the field used for each obser- 
vation. The next column lists the observation dates and the 
last three columns show the filters that were used during each 
observation for each X-ray detector together with the expo- 
sure times after removal of periods of high background. The 
17 XMM-Newton observations gave a total exposure (after re- 
moval of periods of high background) of ~850 ksec for MOS 1 
and MOS2 detectors and ~650 ksec for pn. Some A02 obser- 
vations have an offset of more than ~25 arcmin with respect to 
the other observations. Because of this offset between observa- 
tions, the total exposure time in the centre of the Lockman Hole 
field was reduced significantly, however, the total solid angle 
covered by the observations increased substantially. The solid 
angle of the observation as a function of the effective exposure 
time ' is plotted in Fig.^ 

3. X-ray source list 

We have used the XMM-Newton Science Analysis Software 
(SAS, Gabriel et al. 2004 1 version v5.4, the latest pubHc ver- 
sion of the SAS at the time of study, to analyse the X-ray ob- 
servations. Spurious noise events not created by X-rays were 
filtered from the event files using the XMM-Newton pipeline 

' The solid angle for a value of t was obtained by summing the num- 
ber of pixels in Ml, M2 and pn exposure maps of the total observation 
with a value > t. Because the exposure maps include the energy depen- 
dent mirror vignetting function, the solid angle is given as a function 
of the effective exposure time on each camera 
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rpv/ohs id 




RA 


DEC 




obs. 


date 




Filter'' /GTI'' 








(J2000) 


(J2000) 








EPIC-pn 


EPIC-MOSl 


EPIC-M0S2 


070/0123700101 


PV 


10 52 43.0 


+57 28 


48 


2000- 


04- 


27 


Th/34 


Th/34 


Tck/33 


073 / 01 23700401 


PV 


10 52 43 


+57 28 


48 


2000- 


05 


-02 


Th/14 


Th/14 


Tck/14 


074/0123700901 


PV 


10 52 41.8 


+57 28 


59 


2000- 


05 


-05 


Th/5 


Th/8 


Tck/5 


081/0123701001 


PV 


10 52 41.8 


+57 28 


59 


2000- 


05 


-19 


Th/27 


Th/36 


Tck/28 


345/0022740201 


AOl 


10 5243.0 


+57 28 


48 


2001- 


10 


-27 


M/40 


M/37 


M/24 


349/0022740301 


AOl 


10 5243.0 


+57 28 


48 


2001- 


11 


-04 


M/35 


M/34 


M/31 


522/0147510101 


A02 


10 51 03.4 


+57 27 


50 


2002- 


1 A 


A C 


M/79 


M/81 


M/55 


jz3/014/51Uo01 


AUz 


lU Jl 11 .1 


+57 28 


07 


OA AO 

ZOvZ- 


10 


-17 


M/55 


M/56 


M/4d 


524/0147510901 


A02 


10 52 43.0 


+57 28 


48 


2002- 


10 


-19 


M/55 


M/57 


M/50 


525/0147511001 


A02 


10 52 08.1 


+57 28 


29 


2002- 


10 


-21 


M/78 


M/79 


M/61 


526/0147511101 


A02 


10 53 17.9 


+57 29 


07 


2002- 


10 


-23 


M/45 


M/52 


M/27 


527/0147511201 


A02 


10 53 58.3 


+57 29 


29 


2002- 


10 


-25 


M/30 


M/34 


M/23 


528/0147511301 


A02 


10 54 29.5 


+57 29 


46 


2002- 


10 


-27 


M/28 


M/33 


M/14 


544/0147511601 


A02 


10 52 43.0 


+57 28 


48 


2002- 


11 


-27 


M/104 


M/103 


M/68 


547/0147511701 


A02 


10 52 40.6 


+57 28 


29 


2002- 


12 


-04 


M/98 


M/98 


M/89 


548/0147511801 


A02 


10 52 45.3 


+57 29 


07 


2002- 


12 


-06 


M/86 


M/86 


M/72 



^ Blocking filter: Th: Thin at 40nm Al; M: Medium at 80nm Al; Tck: Thick at 200nm Al 
Exposure time (in ksec) per observation and detector obtained after removal of background flares 




Fig. 1. Solid angle(Q (t)) covered by MOS 1 , MOS2 and pn de- 
tectors as a function of the effective exposure time (after re- 
moval of background flares) in the 0.5-2 keV band. Q. (t) is the 
solid angle (in deg^) in the 0.5-2 keV exposure map of each 
EPIC camera with an exposure >t. 

standard filtering. We have cleaned the event files of the indi- 
vidual observations with periods of high background due to soft 
high background flares. For each observation and detector we 



built light curves (histograms of counts as a function of time). 
The light curves were visually inspected to search for time in- 
tervals affected by high flaring background periods. Thresholds 
in count rate were determined separately for each observation 
and camera. 

Events covering patterns 0-12 for MOS and 0-4 for pn data 
were selected. 

Cleaned event files were used to create images, background 
maps and exposure maps for each detector and for each of 
the XMM-Newton standard energy bands (0.2-0.5, 0.5-2, 2- 
4.5, 4.5-7.5 and 7.5-12 keV). We have run the XMM-Newton 
source detection algorithm, eboxdetect-emldetect, simul- 
taneously in the five energy bands. Our motivation was to reach 
the maximum sensitivity in each individual band in order to 
best detect objects with spectra peaking at different energies 
(as it is the case for AGN with different absorbing column den- 
sities). However, it is important to note that the sensitivity of 
the XMM-Newton X-ray detectors is a strong function of en- 
ergy, with the maximum sensitivity reached between 0.5 and 
4.5 keV. This implies that we will best detect objects with X- 
ray spectra peaking within this interval of energy. 

Due to the large offset between different observations of 
the field (see Sec.|2j, we did not merge the event files, because 
in the merging process important information from the indi- 
vidual observations is lost (e.g., bad columns in the detectors). 
Moreover, it is not possible to create exposure maps or back- 



4 



S. Mateos et al.: X-ray spectra LH 



ground maps from the merged event files. Therefore we ex- 
tracted images, exposure maps and background maps, for each 
individual observation, detector and energy band, and then, we 
combined them to obtain the total observation of the field for 
each X-ray detector and energy band. 

We decided to run the source detection algorithm indepen- 
dently for each detector Because the pn data give the deepest 
observation of the field (the MOS 1 and M0S2 detectors only 
receive about half of the radiation from the X-ray telescopes, 
the other half goes to the Reflection Grating Spectrometers), 
we used the pn source list to build our catalogue of sources. 
However, we cross-correlated the sources detected with the pn 
with the ones detected with each MOS detector. We found that 
only one faint object that was detected with the MOSl and 
MOS2 detectors was not detected with the pn (probably be- 
cause it was very close to a brighter object). We added this 
object to our source list. We have carried out a visual screen- 
ing of the objects to remove spurious detections (e.g. we have 
detections in hot pixels that were still present in the data). The 
final number of objects detected after visual screening in the 
integrated XMM-Newton observation of the Lockman Hole is 
268. 

In this paper we show the results obtained from the spec- 
tral analysis of the 123 sources with the best spectral quality 
(more than 500 MOS-npn background subtracted counts in the 
0.2-12 keV band). Because we were interested in studying the 
X-ray spectral properties of AGN, we excluded from the sam- 
ple the objects identified as clusters of galaxies or stars. The 
results from the analysis of the fainter objects will be presented 
in a forthcoming paper. Fig.|2lshows the distribution of (back- 
ground subtracted) counts for the sources that we have studied. 
In the following we will refer to the sample of 123 brightest 
sources as our list of objects. 

In order to allow comparison with previous surveys con- 
ducted in the soft, 0.5-2 keV, and hard, 2-10 keV, bands, we 
have checked whether our objects were detected in any of these 
bands. We have the likelihoods of detection for the soft band, 
because it was one of the energy bands used for the source de- 
tection. We have calculated the likelihoods in the hard (2-10 
keV) band, combining ^ the ones obtained in the bands 2-4.5 
keV, 4.5-7.5 keV and 7.5-12 keV. Indeed, the combination of 
these values will give us the detection likelihood in band 2-12 
keV. However, the effective area of the X-ray telescopes de- 
creases rapidly at energies above 5 keV, hence we do not ex- 
pect the value of the detection likelihood in the 2-12 keV band 
to differ substantially from the value in the 2-10 keV band. By 
selecting sources with 0.5-2 and 2-10 keV detection likelihoods 
above 10 we found that the number of objects detected in both 
the soft and hard bands was 117 (out of 123). Three objects 
were only detected in the hard band and three only in the soft 
band. These numbers indicate that the source population that 
we are studying is detected in both bands, and therefore our list 
of sources does not differ significantly from what would have 
been detected in a hard or soft band survey at similar fluxes. 




number of counts 



^ Detection likelihoods in the 2-10 keV band were ob- 
tained following the description of the emldetect task, see 
^http://xmm. vilspa.esa.es/external/xmm_sw_cal/sas _frame.shtml^ 



Fig. 2. Distribution of the number of background-subtracted 
counts in the 0.2-12 keV spectra for the sources in our sample. 
We used the counts of the time averaged spectra with better 
signal to noise (MOS or pn) 



4. X-ray spectral products 

An automated procedure has been used to obtain for each in- 
dividual object the spectrum of the total observation, hereafter, 
the time averaged spectrum of the sources. In addition a back- 
ground spectrum and a response calibration matrix have been 
generated by combining individual products from each expo- 
sure. 

First we have extracted the spectra of each object for each 
detector (MOSl, MOS2 and pn) and observation. We used the 
coordinates of the objects, RA and Dec, that we obtained from 
the source detection process, and the SAS task region to de- 
fine the source and background extraction regions. The first 
was defined as a circle with a radius {js) that varied depend- 
ing on the position of the source within the detector. We ex- 
tracted the background for each object in an annulus centred 
on the source position, with inner radius and outer radius 
3xrj. The task region checks the source and background re- 
gions for overlap with neighbouring sources. If overlapping ex- 
ists, then, the size of the source region is reduced until it is 
removed. For the background regions, if neighbouring objects 
fall inside the background region they are masked out. The task 
region also checks that the extraction regions do not extend 
outside the edges of the field of view. The radius of extraction 
of spectra varies from source to source, but typically was ~14- 
20 arcsec. Once the regions were defined, we used the SAS 
task evselect to extract from event files the spectra of each 
object. Calibration matrices (ar£ and rmf) for each spectrum 
were obtained with the SAS tasks arfgen and rmfgen. 

We did not use the spectra from observations where the ob- 
jects were near the borders of the FOV, or near CCD gaps or 
bad columns, because in these cases we found that the spectral 
products, in particular the response matrices were often incor- 
rect. To find and remove these cases we visually checked the 
images of each observation and detector. 
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We have obtained a MOS and pn time averaged spectra for 
each object. As we see in Table [2 different fihers were used 
for the observations. The filters affect in a different way the 
X-ray spectra at low energies. To take into account this effect 
when combining the spectral products of each source, we have 
weighted the data of each individual observation with the ex- 
posure time of the observation. ^ 

The spectra were extracted in the energy range from 0.2 to 
12 keV, where the X-ray detectors are best calibrated. In order 
to use the minimisation during the spectral fitting, we have 
grouped the spectra with a minimum number of 30 counts per 
bin. 



5. Optical identifications 

A large fraction of the X-ray sources that were detected with 
the ROSAT satellite in the Lockman Hole have already been 
identified through optical spectroscopy (Schmidt et al., .1998 
Lehmann et al.. l2001l l. These sources were detected in the 0.5- 
2 keV energy range, and therefore the optical identification is 
expected to be biased against absorbed sources, whose X-ray 
spectra does not peak in the ROSAT energy interval. However, 
the optical identifications go as deep as R~24, hence we do 
not expect the identifications in our sample to be significantly 
affected by the bias against absorbed sources. 

Our XMM-Newton observations cover a larger solid angle 
than the ROSAT observations and find additional sources due 
to the superb high energy response of the XMM-Newton detec- 
tors. Some of the objects that we have analysed fall outside the 
solid angle covered by ROSAT and we do not have optical iden- 
tifications. Other sources were not detected with ROSAT. For 8 
of these newly detected XMM-Newton sources optical spectra 
have been obtained with the LRIS and Deimos instruments at 
the Keck telescopes in 2001, 2003, and 2004 (PI: M. Schmidt 
and P. Henry). The spectroscopic identification of these objects 
and of the entire Deep XMM-Newton Survey in the Lockman 
Hole will be presented in a forthcoming paper (Lehmann et al. 
2005 in preparation). 

To be consistent with the ROSAT identifications we have 
used the same criteria to differentiate between type-1 and type- 
2 AGN as described in detail by Schmidt et al. ( 1998 1. Sources 
were classified depending on the properties of their UV/optical 
emission lines. Objects with UV/optical emission lines with 
FWHM > 1500 km s ' in their optical spectra were classi- 
fied as type-1 AGN. Sources that do not exhibit broad emis- 
sion lines but that show Ne emission lines ([NeV] and/or 
strong [Ne III]) were classified as type-2 AGN. Classification 
of Narrow Line Seyfert 1 galaxies (NLSyl) was only possible 
for bright nearby objects. Hence, we cannot be sure that the 
optical classification criteria used to separate the type-1 and 
type-2 AGN in our sample has excluded all NLSyl from the 
sample of type-2 AGN. However, based on the observed X-ray 
properties of our sources, we do not expect this to affect the 
results of our study. 

' Source and background spectra are obtained adding the counts for 
each channel. The areas used to extract the spectra and the response 
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Fig. 3. Histograms of 0.2-12 keV fluxes for the whole sample 
of sources analysed and for the objects with optical identifica- 
tions. The fluxes were obtained from the sources best fit model 
(See Sec.|9}. 

At the time of this analysis, 74 (~ 60%) of the sources had 
optical spectroscopic identifications. Of these, 46 were classi- 
fied as type-1 AGN and 28 as type-2 AGN. 

In Fig.|3lwe show the distribution in 0.2-12 keV flux of all 
the objects (solid histogram) and of the identified sources (dot- 
dash histogram). We see that the two distributions agree quite 
well, i.e., the identified sources do not tend to have higher X- 
ray fluxes. 

6. X-ray spectral fitting 

We have used the xspec 1 1.3.0 package to fit the X-ray spec- 
tra of our objects. MOS and pn spectra were fitted simultane- 
ously and with the same spectral model, including normalisa- 
tion. At the time of our analysis it was found that there was 
an offset of ~1 arcmin between the optical axis of the three 
EPIC instruments and the values in the calibration files (CCF) 
^. This could be introducing discrepancies in the fluxes mea- 
sured by the MOS and pn due to an incorrect vignetting cor- 
rection. These flux discrepancies could be as high as +14%. 
However, we did not find a significant improvement in the qual- 
ity of our fits when different normalisations were used to fit 
MOS and pn spectra. For the objects where we found an off- 
set between MOS and pn normalisations, this offset was much 
higher than the expected flux discrepancies explained before. 
We interpreted this effect as a change in the flux of the source 
during the observations (note that the time averaged spectra of 
MOS and pn in the majority of the cases were not necessary 
built with the same set of observations, because the FOV and 
the positions of the gaps are different for MOS and pn, and in 
different observations). 

matrices were weighted with the exposure times of each individual 
observation. 

* Fluxes were obtained from the objects best fit model (See Sec.|9} 
^ http://xmm.vilspa.esa.es/docs/documents/CAL-SRN-0156-l- 
S.ps.gz 
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Fig. 4. Dependence of F with soft (0.5-2 keV) and hard (2-10 keV) flux when ah spectra are fitted with a single power law model. 
In plots (a) and (b) we show the dependence of the weighted (with the error of each individual value) F with the flux. The bins in 
flux were defined in order to have the same number of sources per bin. In plots (c) and (d) we show the values of F obtained for 
each individual source. The dash-dot lines in these plots show, for an exposure time of 280 ksec, the limit in flux for detection 
for an object as a function of F (see Sec. l7.1l for details). Errors bars in (c) and (d) correspond to 90% confidence. 



In order to compare the results of our analysis with other 
studies of data with lower signal to noise, we have fitted the 
spectra of our sources with a single power law model (hereafter 
SPL). This has allowed us to study in more detail the origin 
of the hardening of the average spectra of AGN with the soft 
X-ray flux (see e.g. Mateos et al 2005 1, and investigate if the 
same eff'ect is also present when 2-10 keV fluxes are used. To 
study the eff'ect of absorption in the results obtained with the 
SPL we fitted all the spectra with an absorbed power law model 
(hereafter APL model). 

Using these two models it is possible to obtain very use- 
ful results on the average spectral properties of our sources, 
the broad band continuum shape and the X-ray absorption. 
However, our major goal is to study in detail the 0.2-12 keV 
X-ray emission of each individual source. Hence, for each ob- 
ject we have obtained its best fit model. The quality of our data 
has allowed us to search for soft excess at low energies, the Fe 
Ka emission fine complex, and reflection components. 



We have used the F-test to measure the significance of de- 
tection for each spectral component. We have selected a con- 
fidence level threshold of 95% to accept an additional spectral 
component as being real. The criteria that we used to select X- 
ray absorbed sources is to have an F-test significance > 95%. 
This is different from some definitions found in the literature, 
because we did not impose a lower threshold in the detected 
values of Nh. For example, Ueda et al. ( I2003> defined as X- 
ray absorbed the sources with absorption column density at the 
source redshift > 10^^ cm"^. However, it is important to note 
that with our criteria all sources selected as X-ray absorbed had 
values of Nh > 21 cm"^. 



7. Single power law fitting (SPL) 

We have used a single power law model to fit the 0.2-12 keV 
emission of all the objects. In this model, the free parameters 
are the normalisation (the same for MOS and pn spectra) and 
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the slope of the broad band continuum, F ^. The power law is 
absorbed with a fixed column density of 5.7 x 10'^ cm"^ to in- 
clude the effect of the absorption by our Galaxy in the direction 
of the Lockman Hole field. 

7.1. Dependence of V with the X-ray flux 

The results of the fits have been used to study the dependence 
of r with the X-ray fluxes obtained from the SPL model. The 
results are plotted in Fig. |2 In plots (a) and (b) we show the 
dependence of (F) with the 0.5-2 keV and 2-10 keV fluxes. 
The bin sizes were defined in order to have the same number of 
objects per bin, and the average values were obtained weighting 
with the eiTors of each individual value. When a single power 
law model is used, we see that the average continuum shape 
becomes harder with decreasing 0.5-2 keV flux. However, it is 
interesting to note that we do not see any dependence of (F) 
with the 2-10 keV flux down to ~ 3 x 10"'^ erg cm"^ s"'. 

The values of (F) were calculated with the standard formula 
for the weighted mean, 

(F> = ^ P,- X F; 

where the weight, P,-, of each individual best fit value, F,, is a 
function of the error in the parameter obtained from the fit, o",, 
i.e. 

p. - ; ! 

' E(i/^i=) 

To calculate the uncertainty in (F) we have used the error on 
the weighted mean (Bevington et al. 1992j, 

that includes the measurement errors, cr,, and the dispersion of 
each F, from the estimated value (F). Using these expressions 
we found that our objects have (F)= 1 .79+0.03 when their spec- 
tra are fitted with a SPL (the value is 1 .60+0.05 if we use the 
unweighted mean). 

In order to understand better the origin of the hardening of 
(F) with the 0.5-2 keV flux, and why we do not see the same ef- 
fect using 2-10 keV fluxes, we have plotted in Fig.|3(figures (c) 
and (d)) the values of F that we obtained for each individual ob- 
ject. Thanks to the good quality of our data we can see that (F) 
becomes harder because at faint 0.5-2 keV fluxes a population 
of faint sources is revealed with very hard (< 1) spectral slopes. 
We also see that the number of faint hard objects becomes more 
important as we go to fainter fluxes. In the 2-10 keV band we 
do not see the hardening of (F) because these hard objects are 
detected at all 2-10 keV fluxes. Moreover, their number seems 
not to vary with the 2-10 keV flux. 

We have studied whether our criteria for selection of ob- 
jects (i.e. MOS+pn background subtracted counts above 500) 
could be introducing any bias in our results. In particular we 
have studied whether, for a given flux, we are favouring ob- 
jects with a given spectral slope. To study this, we have caiTied 

^ We use the power law photon number index, F. Its relation with 
the energy index is a = F-l 



out simulations. We first defined a grid of points in F and S 
(first using 0.5-2 keV flux and later with 2-10 keV flux), cov- 
ering the same range of F-S values as our sources. Using a pair 
of on-axis response matrices, arf and rmf, and typical back- 
ground spectra selected from one of our objects, we have simu- 
lated a spectrum on each grid point. With the simulated spectra 
we have calculated the minimum exposure time that is needed 
to reach the threshold in number of counts that we have used 
to select our sources (i.e. 500 MOS+pn background subtracted 
counts). With these simulations we were not interested in quan- 
tifying the limits of detection as a function of S and F, but to 
study the biases in our sample, and whether they affect our re- 
sults. Therefore we only need to do one simulation on each grid 
point and then we just have to search for the points in the F-S 
grid with the same value of the exposure time. 

A constant exposure time line to get 500 counts is repre- 
sented with the dot-dash lines in plots (c) and (d) on Fig. 0] 
for an exposure time of 280 ksec. We see that in the soft 
band, for a given flux, we have the same efficiency of detec- 
tion for different values of F down to 1.5. At fluxes above 
~ 10"'^ erg cm"^ s"' this bias is not affecting the observed 
hardening of F. Only the bin at the faintest 0.5-2 keV fluxes 
(plot (a)) could be affected by this bias. Our simulations show 
that the hardening in (F) is an intrinsic property of our objects. 
As we said in Sec. |3l our objects were all detected in the soft 
band, therefore this effect is a property of the 0.5-2 keV popu- 
lation of objects. The objects responsible for the hardening of 
(F) can be more absorbed sources or sources having intrinsi- 
cally harder spectra. 

In the 2-10 keV band we obtained different results from 
the simulations. At the faintest fluxes we most easily detect 
objects with soft spectra. This is an expected result because 
the effective area of the X-ray detectors in XMM-Newton de- 
creases rapidly at energies >5 keV and therefore it is more dif- 
ficult to detect faint objects with flat spectral slopes. However, 
down to the flux level where we start to lose faint hard objects, 
~ 6 X 10"'^ erg cm"^ s"', we see that there is no dependence of 
F with hard flux because hard objects are detected at all flux 
levels. 



8. Absorbed power law fitting (APL) 

In order to study the nature of the population of faint hard 
sources responsible for the hardening of {F) with the 0.5-2 keV 
flux we have fitted the spectra of all the objects with an ab- 
sorbed power law model. The free parameters of this model 
are the normalisation, the spectral slope of the power law com- 
ponent, and the intrinsic (rest-frame) absorption of the objects 
that we have measured in the observer's frame (N*'^). Again, 
we also included the effect of the Galactic absorption in the 
direction of the Lockman Hole. Using the APL model we ob- 
tained a value for the weighted mean of (F) - 1.87 + 0.04 (the 
value being (F) = 1.95 + 0.08 using the arithmetic mean) for 
the objects where absorption was detected (F-test > 95%), and 
<F) = 1.95 + 0.03 (the value being <F) = 1.82 + 0.04 using 
the arithmetic mean) for the objects where we did not detect 
absorption (F-test < 95%). 
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Fig. 5. Dependence of (F) with 0.5-2 and 2-10 keV flux for absorbed (F-test > 95%) and unabsorbed (F-test < 95%) sources. For 
each source we used F and S from its best fit model (single power law or absorbed power law). See Sec.|Slfor details. 



The dependence of (F) on 0.5-2 keV and 2-10 keV fluxes 
that we see fitting the spectra of our objects with the APL 
model is shown in Fig.|5l For the objects where absorption was 
detected we used the spectral parameters (F and observed S) 
obtained with the APL model. For the unabsorbed sources we 
used the parameters from the SPL model. We see that absorp- 
tion can account for most of the hardening of (F) with the soft 
X-ray flux. We also see in Fig.|5]that the same dependence of 
(F) with the X-ray flux is obtained for absorbed and unabsorbed 
objects. 

Mateos et al. ('20051 found that for their serendipitous X-ray 
sources (with much lower spectral quality), the average spec- 
trum of unabsorbed (F-test significance <95%) sources signif- 
icantly hardens at faint 0.5-2 keV fluxes. They concluded that 
undetected absorption was responsible for the observed effect 
in unabsorbed sources. The magnitude of undetected absorp- 
tion cannot be very significant in our Lockman Hole sources as 
we do not see a clear hardening of (F) for unabsorbed sources. 
Using F and S from the APL model for unabsorbed sources 
we obtain the same result. Hence, if we still have sources with 
undetected absorption, their absorbing column densities cannot 
be very high. Moreover, we do not find evidence for the exis- 
tence of a population of faint sources with intrinsically harder 
spectral slopes. 

The results for the hard fluxes do not vary significantly from 
what we obtained with the SPL fits. This is exactly what we 
would expect if absorption produces the hardening of F, be- 
cause the 2-10 keV fluxes are less affected by absorption. 

In Fig.|6lwe plot the values of absorption (observer's frame) 
that we obtained with the APL model, as a function of the 

0. 5-2 and 2-10 keV fluxes. There is an obvious correlation be- 
tween the absorption and the observed soft flux. We see that 
the distribution of absorbing column densities does not seem 
to vary with the hard band fluxes. We have studied the depen- 
dence of N*** with the 0.5-2 ke V flux using de-absorbed fluxes, 

1. e. the 0.5-2 keV fluxes corrected for the effect of absorption. 
We found that when de-absorbed fluxes are used, N?}''* does not 



vary with X-ray flux, i.e., fainter and/or more distant objects do 
not appear to be more absorbed. 

We have studied whether the fraction of X-ray absorbed 
objects depends on the flux after correcting for the effect of 
absorption (the Nh columns measured in our sample of sources 
are not high enough as to affect significantly the measured 2-10 
keV fluxes and hence the correction of 2-10 keV fluxes for the 
effect of absorption is not significant). The results are plotted 
in Fig0 For comparison we have plotted the results that we 
obtain when absorbed fluxes are used (circles). We do not see 
significant differences using absorbed or de-absorbed 2-10 keV 
fluxes, because as explained before, these are not significantly 
affected by the absorption measured in our sources. However 
important differences are seen when 0.5-2 keV fluxes are used. 
If we do not correct for the effect of absorption in the 0.5-2 keV 
flux, we see an increase in the fraction of absorbed objects at 
fainter fluxes. However, if de-absorbed fluxes are used instead, 
the fraction of absorbed objects does not vary with the X-ray 
flux and we obtain the same result as in the 2-10 keV band. 

Note that the fraction of absorbed sources at the faintest 
de-absorbed 0.5-2 keV fluxes (~ 20%) is significantly lower 
than the values found at brighter fluxes. The absorbed sources 
that should contribute to the bin at the faintest fluxes have an 
observed flux below the threshold applied to our objects and 
therefore are not included in our sample (remember that to se- 
lect our sources we used 0.2-12 keV counts, i.e. ~fluxes with- 
out correction for absorption). Another effect that could also 
contribute to this result is that at the faintest fluxes we may 
have some sources with undetected absorption. We do not ex- 
pect this effect to be important, because as we see in Fig.|5l{F) 
for unabsorbed sources does not seem to become significantly 
harder at the faintest 0.5-2 keV fluxes. 

9. Best fit model 

Up to now we have shown the results from spectral fits where 
only the X-ray continuum and the intrinsic absorption were 
modelled. However, there are other spectral components that 
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Fig. 6. Distribution of N"*"* for absorbed (F-test > 95%; squares) and unabsorbed (F-test < 95%; triangles) objects. Error bars 
correspond to 90% confidence. 
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Fig. 7. Fraction of absorbed (F-test>95%) objects vs 0.5-2 (left) and 2-10 (right) keV flux. For each source we used the fluxes 
from its best fit model (single power law or absorbed power law). Circles show the results that we obtained using observed fluxes, 
while the triangles show the results using fluxes corrected for the effect of absorption. See Sec.|8lfor details. 



Table 2. Results from the X-ray spectral analysis. 



Model" 


Total 


type-1 AGN 


type-2 AGN 


Not id.'' 


SPL 


65 


35 


5 


25 


SPL + SE 


4 


4 








APL 


39 


6 


16 


17 


APL + SE 


5 





4 


1 


CAPL 


9 


1 


3 


5 


2SPL 


1 








1 


Total 


123 


46 


28 


49 



" Best fit model: SPL: single power law; APL: Absorbed power law; 
SE: Soft excess; CAPL: partial covering; 2SPL: two power laws 
(see Sec.|21for details). 
Objects without optical identifications. 



can also contribute significantly to the emission in the 0.2-12 
keV energy band. There is evidence for them in the results that 



we have shown previously. For example we see in Fig. |5l that 
(F) seems not to vary with the X-ray flux when absorption is in- 
cluded in the fitting model. However we still have a clear scatter 
in the points which cannot be explained if F does not vary sig- 
nificantly with flux for the objects in our sample as our results 
appear to show. However, we would expect this scatter of (F) if 
other spectral components are present in the data (e.g. soft ex- 
cess emission) and they are not properly modelled. Besides the 
soft excess, other spectral components that can contribute to the 
0.2-12 keV emission are ionised absorption, the Fe Ka emis- 
sion complex and the Compton reflection hump that should ap- 
pear at high X-ray energies. 

We have studied in detail the MOS and pn time averaged 
spectra of each individual object. To model the soft excess 
emission we have used a black-body model. This component 
adds two free parameters to the fit, the temperature (in keV) and 
normalisation of the black-body. For some objects we could not 
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get a good fit of the detected soft excess with a black-body. In 
all these cases we obtained a good fit using a partial covering 
model (i.e., only part of the X-ray emission from the inner most 
region of the AGN is absorbed) to fit the signatures of absorp- 
tion and soft excess emission. This model introduces one new 
parameter to the fit with respect to the APL model, the covering 
fraction of the absorber (between and 1). The only emission 
line that we expect to detect with the quality of our spectra is 
the Fe Ka complex at 6.4 keV (rest frame energy for neutral 
iron). To search for this component we have used a Gaussian 
line profile, that allows us to calculate the centroid (in most 
cases we fixed the centroid to 6.4 keV), width and normali- 
sation of the line. Absorption signatures found in some spectra 
were modelled with an absorption edge (zedge in xspec). This 
model introduced two free spectral parameters, the threshold 
energy and the absorption depth at the threshold energy. The 
most prominent signature from reflection in AGN is a change 
in the slope of the X-ray continuum at energies above 10 keV 
(rest frame). This component is known as the Compton reflec- 
tion hump. We do not expect our objects to be bright enough as 
to detect with high significance reflection signatures given the 
limited bandpass of XMM-Newton. However, we have searched 
for this component in all the spectra adding a second power law 
to the fits at high energies. 

The best fit model for each source is the one that gave a 
significant improvement in fit over the previous one in the se- 
quence SPL-APL (with z in the case of identified sources)- 
APL-Hsoft excess. The improvement in the fit was measured 
by the usual F-test, taking into account the improvement in the 

value and the number of new parameters introduced. For 
some objects with soft excess emission, the detection of ab- 
sorption was only significant after modelling of this spectral 
component. In particular, when the classification (type-l/type- 
2 AGN) of the identified sources is used, the parameters from 
flie APL model with intrinsic absorption are always used. 

Our results are listed in Table |2 The best fit spectral pa- 
rameters obtained for each object are shown in Tabled In the 
present paper we identify our objects with the numbers that 
will be used in a forthcoming catalogue paper of the Lockman 
Hole (H. Brunner et al. 2005, in preparation). The results were 
obtained fitting MOS and pn spectra simultaneously. If it was 
not required by the data we used the same MOS and pn model 
normalisations. For the sources for which we used different 
MOS and pn normalisations we show fluxes and luminosities 
obtained with the spectrum with best quality. 

For 65 objects (~53%) the best fit model was a single power 
law. Among these objects there were 35 type-1 AGN and 5 
type-2 AGN. We detected absorption in 53 objects out of 123 
(43%), including 7 type-1 AGN and 23 type-2 AGN. In Table|3 
we list the fraction of absorbed objects in our sample and in the 
samples of type-1 and type-2 AGN. The values were corrected 
for spurious detections ^ following the method described in 
Mateos et al. (2005 1. We used the method described in Stevens 



et al. J2005> and Mateos et al. J2005> to compare the fraction of 
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Fig. 8. Dependence (F) with 0.5-2 keV flux. For each source 
we used the spectral parameters obtained from its best fit model 
(see Sec.|9lfor details). 



absorbed objects in type-1 and type-2 AGN. We find that these 
fractions are different (the number of absorbed objects in type- 
2 AGN being larger than in type-1 AGN) with a significance of 
more than 99.99%. 



Table 3. Results of detection of X-ray absorption. 





N.O, 


Nabs 


r 


■> lint 


All sources 


123 


53 


0.38 


>0.27 


type-1 AGN 


46 


7 


0.10 


<0.29 


type-2 AGN 


28 


23 


0.77 


>0.51 



' We accepted the spectral signatures as being real if the signifi- 
cance of detection from the F-test was >95%. Hence 5% of detections 
are expected to be spurious. 



" Fraction of absorbed objects taking into account the expected frac- 
tion of spurious detections (see Sec|5) 
^ 3cr limits in the fraction of absorbed objects 



9.1. Broad band continuum 

Using the best fit spectral slopes we obtained a weighted (F) 
for the objects in our sample of 1.87 + 0.02 (1.86 + 0.02 if the 
arithmetic mean is used). In Fig.|Slwe show the dependence of 
(F) with the 0.5-2 keV flux using the best fit parameters (F and 
S) for each object (the same result is obtained using 2-10 keV 
fluxes). We see that when we take into account during fitting 
all the spectral components, the scatter in (F) (see Fig. |5j is 
reduced significantly, and most of the points are consistent with 
the obtained average value of F. 

The mean spectral slope was found to be 1.89 ±0.03 (1.88 + 
0.03 with the arithmetic mean) for type-1 AGN and 1.71 +0.03 
( 1 .82+0.06 with the arithmetic mean) for type-2 AGN. In Fig.|9l 
we compare (F) for type-1 and type-2 AGN as a function of the 
X-ray fluxes. For type-1 AGN, where only 7 out of 46 are X-ray 
absorbed, we obtain the same results using 0.5-2 and 2-10 keV 
fluxes, i.e., no dependence of (F) with the X-ray flux. However, 
X-ray absorption is important in type-2 AGN (>50% of type-2 
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Fig. 9. Dependence of (F) with 0.5-2 and 2-10 keV flux for type-1 and type-2 AGN. For each source we used F and S from its 
best fit model (single power law or absorbed power law). See Sec. l9.1l for details. 



AGN being X-ray absorbed). The ratio of type-2 AGN/type-1 
AGN increases as we go to fainter 0.5-2 keV absorbed (i.e. not 
absorption corrected) fluxes, while it remains constant with 2- 
10 keV flux. We explained in detail in Sec.|6lthat this is due to 
the existing correlation between Nh and the soft absorbed flux, 
i.e. most absorbed objects have the faintest fluxes in the 0.5-2 
keV band. We see in Fig.^that there is a clear dispersion in (F) 
for type-1 and type-2 AGN, however it seems that type-2 AGN 
tend to have lower (F) than type-1 AGN at the fluxes covered 
by our sample. Excluding from the sample of type-2 AGN the 
three sources that we found with no detected X-ray absorption 
and spectral slope significantly lower than the average value 
for type-1 AGN (see Sec. II Hand Table^}, we still see the same 
dependence of (F) with X-ray flux for type-2 AGN. 

We have followed the procedure described in Nandra & 
Pounds (T994'l and Maccacaro et al. (1988) to estimate the 
intrinsic dispersion of the photon index in type-1 and type-2 
AGN, and to investigate whether after allowing for intrinsic 
dispersion in F, we still find type-2 AGN to be on average flatter 
than type-1 AGN. In this method it is assumed that the disper- 
sion in F values can be described well with a Gaussian function 
of mean (F) and dispersion cr^p). The results of this analysis are 
listed in Table |3 where we have the values obtained using the 
weighted mean for comparison. We have found that there is an 
intrinsic dispersion in F of ~ 0.2 in type-1 AGN and type-2 
AGN, and that the value of the dispersion is similar in both 
samples of objects. It is interesting to note that the results ob- 
tained with the weighted mean and the Maximum Likelihood 
method are consistent within each other. 

In Fig.[lOlwe show the contours in (F)-cr(p) space for a Ax^ 
of 2.3, 6.17 and 11.8 that correspond to 1, 2 and 3cr for two 
parameters. The significance of type-2 AGN being on average 
flatter than type-1 AGN is only at 1 .62cr (using the values of F 
obtained with the ML method ^). 



* Note that the significance of F being different for type- 1 AGN and 
type-2 AGN is ~ 4cr if the values obtained with the weighted mean 
are used and no intrinsic dispersion in F is considered 



However, it is important to note that if the signatures of 
absorption in the X-ray spectra are not very significant, the de- 
tected values of Nh will tend to be lower than the real ones 
(see Mateos et al. I2005> and then, the fitted F wiU be flatter 
We expect this effect to be more important for type-2 AGN 
where we have more sources with absorption. The small differ- 
ence in (F) for type-1 and type-2 AGN might be due to this ef- 
fect. Therefore with the current data we cannot reach any strong 
conclusion. 

9.2. X-ray absorption 

We have detected X-ray absorption in ~ 37% of the sources in 
our sample. Absorption was found in ~ 10% of type-1 AGN 
and ~ 77% of type-2 AGN. We first checked that the measured 
F and Nh were not correlated, and therefore that we have ob- 
tained reliable parameters, especially the column density for 
each individual object. The results are plotted in Fig. for 
sources with known redshifts, where we do not see any evident 
correlation between the two spectral parameters. In the objects 
with large F, the values including the error bars are in all cases 
consistent with a value of F ~2. 

Table 4. Mean spectral photon index of type-1 and type-2 AGN 
obtained with the weighted and arithmetic means and with the 
Maximum Likelihood analysis. The spectral slopes from the 
sources' best fit model were used. 



Maximum Weighted Arithmetic 
Likelihood Mean Mean 
Sample (F) tr 

Whole sample 1.92j;'j'^ 0.28j;'j'^ 1.87 ±0.02 1.86 ±0.02 

type-1 AGN 1.89;!-!J^ 0.20»-«^ 1.89 ±0.03 1.88 ±0.03 

type -2 AGN 1.7255-^^ 0.20°-'° 1.71 ±0.03 1.82 ±0.06 
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Fig. 10. Contour diagrams for the value of the average spectral 
slope and intrinsic dispersion of our samples of type-1 AGN 
(solid lines) and type-2 AGN (dashed lines) obtained from the 
Maximum Likelihood analysis (see Sec. |9jl. The contours are 
defined as b^x"-'^-^^ 6.17 and 11.8 corresponding to standard 
1, 2 and 3cr confidence regions for two parameters. 




Fig. 11. r vs. NJ^"^ (rest-frame) for the type-1 and type-2 AGN 
with detected absorption. Note that all type-1 AGN have col- 
umn densities between 10^' - 10^^ cm"^ while type-2 AGN 
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Fig. 12. Distributions of intrinsic (rest frame) absorption in 
type-1 and type-2 AGN obtained from the best fit model for 
sources. 



A type-1 AGN unabsorbed 

A type-1 AGN X-ray obsorbec 

□ type-2 AGN unabsorbed 

■ type-2 AGN obsorbed 



■ A 



2 

redshift 



Fig. 13. 2-10 keV luminosity (corrected for absorption) vs. red- 
shift for type-1 and type-2 AGN. Absorbed sources are the ob- 
jects where we found an F-test significance of improvement of 
the fits > 95%. 



have a much wider distribution of Nj^"". 



Note in Fig.^] that the NJ^"^ distributions in type-1 and 
type-2 AGN seem to be diff'erent. The measured column den- 
sities in absorbed type-1 AGN are between lO^^ - lO^^ 

cm , 

while type-2 AGN have a much wider distribution of values, 
many objects having Nj^" > 10^^ cm"-. We show the distribu- 
tions of Nj!|" in type-1 and type-2 AGN in Fig. [121 The distribu- 
tions appear to be different, with type-2 AGN being in general 
more absorbed than type-1 AGN. Using the KS test to compare 
the two distributions we obtained a probability of them being 
different of > 92%. We will have to wait for the analysis of the 
faint sample of objects before reaching a stronger conclusion. 

In terms of the unified model of AGN, X-ray absorption and 
optical obscuration should be correlated. There is observational 



evidence that this does not hold for all AGN, although only 
a few cases are based on a proper spectral analysis. We can 
confirm that these discrepancies exist in a significant fraction 
of our objects. Mateos et al. (jJOOS ) found some indications that 
the AGN/host galaxy contrast effect might explain why broad 
optical lines are not observed in unabsorbed type-2 AGN. In 
Fig. [0]we show the redshift distribution of the AGN in our 
sample. Because our sample of AGN is much smaller we can 
not perform a similar test. 



' In the redshift interval where the distributions of type-1 and type- 
2 AGN overlapped, they found that unabsorbed type-2 AGN were less 
luminous than type-1 AGN 
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Table 5. Properties of the soft excess emission in type-1 and type-2 AGN that was modelled with a black body. 
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ID 


Class 


redshift 


(Lbb/Lpo) 


kT 


Lbb 


L 










(0.5 - 2keV) 


(eV) 


(0.5 - 2 keV) 


(2 - lOkeV) 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


90 


type 


- 1 AGN 


0.467 


0.363 


105!i^ 


42.64 


43.08 


148 


type 


- 1 AGN 


1.113 


0.479 


78+f 


44.44 


44.74 


270 


type 


- 1 AGN 


1.568 


0.426 


ior\° 


43.80 


44.12 


342 


type 


- 1 AGN 


0.586 


0.081 




42.79 


43.78 


259 


type 


-2 AGN 


0.792 


0.063 


9Q1+149 


42.11 


43.90 


290 


type 


-2 AGN 


0.204 


0.035 


188^;^ 


41.32 


42.83 


424 


type 


-2 AGN 


0.707 


0.060 


474+26 

^' -153 


42.25 


43.74 


511 


type 


-2 AGN 


0.704 


0.014 


83!38 


42.12 


43.58 



Columns are as follows: (1) Source X-ray identification number; (2) object class based on optical spectroscopy; (3) redshift; (4) ratio of soft 
excess to power law 0.5-2 keV luminosities (this ratio is frequently used to measure the strength of the soft excess emission); (5) temperature 
of the soft excess (using a black-body model); (6) logarithm of the 0.5-2 keV luminosity of the soft excess component; (7) logarithm of the 
2-10 keV luminosity of the power law component (for absorbed sources the luminosity was absorption corrected). 



9.3. Soft excess 

We detected soft excess emission (>95% confidence limit from 
an F-test) in 18 (15%) objects. The number of MOS+pn counts 
in the soft excess component vary from 100 to 1000 except for 
one source where the soft excess component has ~3000 counts. 
Within the sources classified as AGN we found soft excess in 5 
(11%) type-1 AGN and 7 (25%) type-2 AGN. The significance 
of the fractions of type-1 and type-2 AGN with detected soft 
excess emission being different is 97%. 

Although our results suggest that soft excess emission 
might be more common in type-2 AGN, it is important to note 
that our samples of type-1 and type-2 AGN have different red- 
shift distributions, and for the highest redshift sources (all type- 
1 AGN) we expect most of the signatures of soft excess emis- 
sion to be redshifted outside the observed energy interval, mak- 
ing the detection of soft excess more difficult. Therefore, we 
have repeated the comparison using only sources in the red- 
shift interval where we detected soft excess (z< 1.568). In this 
case, the significance of the fractions of type-1 and type-2 AGN 
with detected soft excess emission being different is reduced to 
87%. Hence, with our data we cannot confirm that soft excess 
emission is more common in type-2 AGN than in type-1 AGN. 

In 9 sources (4 type-1 AGN, 4 type-2 AGN and 1 unidenti- 
fied object) we fitted the soft excess emission with a black body 
model (a Raymond Smith model gave an equally good fit). The 
properties of the soft excess, i.e. temperature, 0.5-2 keV lu- 
minosity and strength, for the identified sources are listed in 
Table |5] We see that the measured black-body properties do 
not depend on the 2-10 keV X-ray luminosity of the objects. 
For the unidentified object with detected soft excess fitted with 
a black body we found an observed black body temperature of 

The average temperature of the black body was found to be 
0.09±0.01 keV for type-1 AGN and 0.26+0.08 keV for type-2 
AGN. The average 0.5-2 keV luminosities of the black body 
(in log units) were 43.42 + 0.43 erg s"' for type-1 AGN and 
44.1 1 + 0.44 erg s"' for type-2 AGN. The 0.5-2 keV luminosi- 
ties of the soft excess component in type-1 and type-2 AGN 



were not found to be significantly different (a KS test of the 
luminosity distributions gave a significance of them being dif- 
ferent of only 90%). However the measured temperatures of 
the soft excess were found to be higher in type-2 AGN than in 
type-1 AGN. This could be because the soft excess in type-2 
AGN might contain a fraction of scattered radiation. We also 
see that in most sources the temperatures of the black body 
are well above 60 eV and hence it is difficult to explain the 
origin of the soft excess emission in these sources as thermal 
emission from the accretion disc. Componization of cool pho- 
tons in a cloud of hot electrons surrounding the accretion disc 
might be an alternative explanation. 

In 9 other objects (1 type-1 AGN, 3 type-2 AGN and 5 
unidentified objects) the black body could not fit the signa- 
tures of the soft excess. An alternative method for modelling 
the curvature at soft energies is a scattering or partial cover- 
ing model (pcf abs in xspec). The model consists of the sum 
of two power law components having the same spectral index, 
but affected by different absorption (quantified with the cov- 
ering fraction parameter). This model improved significantly 
the quality of the fits, and provided a good fit of the soft ex- 
cess emission in all sources. The average covering fraction that 
we obtained was 0.82+0.06 (the maximum and minimum val- 
ues being 0.98 and 0.50). This value implies that the scattering 
fraction in these sources is rather large (18 + 6%). 

9.4. Reprocessed components 

We have searched for a flattening of the continuum at high 
energies (i.e. Compton reflection) adding a second power law 
to the model. We have found signatures of spectral hardening 
at high energies in only one object, source 108, which is still 
unidentified. The unfolded MOS and pn spectra of this source 
are shown in Fig.^^ We first fitted the X-ray spectrum of this 
object with a single power law giving F ~1.5 but the fit was 
poor, with ax^ of 127 for 67 degrees of freedom. We found that 

The hottest thermal emission expected from an accretion disc sur- 
rounding a 10* Mq black hole. 
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Table 6. Parameters of the Gaussian line in the six identified sources where signatures of line emission where detected with an 
F-test significance > 95%. 



ID 


Class 


redshift 


L2-10 


Line 


<T 


Equivalent 


F - test 










erg s"' 


energy (keV) 


(keV) 


width (eV) prob. (%) 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


270 
21 


type 
type 


- 1 AGN 
-2 AGN 


1.568 
0.498 


44.12 
43.28 


5 90+" '^ 


00+" '*- 

"•""-0.00 

"•-^"-0.45 


452 
1462 


99 
95 


172 
290 


type 
type 


-2 AGN 
-2 AGN 


1.170 
0.204 


44.00 
42.83 


6 40+187 

"■-'^-0.61 


09+^ '^ 

"•"=^-0.00 

00+" " 
"■""-0.00 


360 
283 


99 
96 


326 
407 


type 
type 


-2 AGN 
-2 AGN 


0.780 
0.990 


44.08 
44.20 


6 40+" '" 


02+^ '- 
"■"^-0.00 

40+"'" 

"•^"-0.18 


224 
653 


99 
> 99.99 



Columns are as follows: (1) Source X-ray identification number; (2) optical class based on optical spectroscopy; (3) redshift; (4) 2-10 keV lumi- 
nosity; (5) and (6) rest-frame centroid energy and width of the emission line; (7) rest-frame equivalent width of the line; (8) F-test significance 
for the detection of the line. 




1 2 
channel energy (keV) 

Fig. 14. Unfolded MOS and pn spectrum of the only source in 
the sample where we detected a hardening in the X-ray con- 
tinuum emission at high energies (Compton reflection hump). 
Because this object is still unidentified we fitted the hard spec- 
tral component with a second power law. The F-test signifi- 
cance of this component was found to be 99.98%. 



there was a clear excess emission at high energies. We then fit- 
ted the spectrum with two power laws, and the significantly 
decreased to 98 for 65 degrees of freedom. The F-test signifi- 
cance of improvement of the fit with the new component was 
99.98%. In this case we obtained a value of F of 1.83 ± 0.17 
for the continuum emission flattening out to F = -2.56^Q jg at 
high energies. We did not find evidence for X-ray absorption 
or emission lines in the spectrum of this source. We will have 
to wait until we have the optical identification of this source 
before saying more about the X-ray emission of this object. 

Another signature of reprocessing that has been found in 
many spectra of AGN is an emission line around 6.4 keV. 
This is interpreted as Fe Ka fluorescence from cold matter 
(Pounds et al., I1989I1990I Nandra et al., [79911 Nandra & 
Pounds 119941 and might originate from the reprocessing of 
hard X-ray photons in the accretion disc (Pounds et al. .l990j . 



We have searched for this component in our sources using a 
Gaussian model. Most time averaged spectra of our sources 
do not have enough signal to noise as to detect Fe Ka line 
emission. However, we have been able to detect signatures of 
line emission with an F-test significance > 95% in the MOS 
and pn spectra of 8 objects (1 type-1 AGN, 5 type-2 AGN 
and 2 unidentified sources). In Fig ^] we show the MOS-i-pn 
unfolded time averaged spectra of these sources. The parame- 
ters of the Gaussian line for each identified source are listed in 
Tableil 

Sources 21 and 407: In these sources we found a signif- 
icant width in the line profile, which might be indicating that 
the line was formed in the inner parts of the accretion disc, 
and hence it should have a relativistic profile (with a red wing 
component due to gravitational redshift). When fitted with a 
Gaussian model, we would expect the line centroid to be found 
at an energy slightly below 6.4 keV. While in source 407 the 
line centroid is consistent with being neutral iron, in source 21 
it was ~5.9 keV (although consistent with being neutral iron 
within the error bars). 

Sources 270 and 290: In these sources we detected a nar- 
row Gaussian line and line centroids lower (but consistent 
within the error bars) than the value for neutral iron. 

Sources 172 and 326: In these sources we also found a 
significant line width, although in both cases it was consistent 
with zero at 90% confidence. In source 172 the line centroid 
was consistent with being neutral iron, but in source 326 it was 
significantly higher (even within the error bars). In this source 
the line might be originating in an ionised accretion disc. 

In all the spectra where we detected the line we did not 
have enough signal to noise in the data in order to use a more 
physical model to fit the profile of the line (xspec models laor 
for a Kerr black hole or diskline for a Schwarzchild black 
hole). The rest frame equivalent width (EW) of the line in the 
type-1 AGN where we detected this component was found to 
be ~452 eV. In most type-2 AGN the measured values were 
between 200 and 600 eV. However there is one source, 21 for 
which we found a rest frame EW of ~1400, substantially higher 
than in the other type-2 AGN. It is important to note that in this 
source the F-test significance of detection of line was the lowest 
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Fig. 15. Unfolded MOS and pn time averaged spectra of the 6 AGN where we detected signatures of emission line at high energies 
(F-test significance >95%). In the X-ray spectrum of source 172 (type-2 AGN) we also found an absorption edge at an energy of 
~ 7.56!!^g keV with absorption depth t - l-4g ? (F-test significance of detection was 99%)). 



among all sources (95%) and hence the measured value has the 
highest uncertainty. 

Using the same sample of objects, Streblyanska et 
al. (2005) found a clear relativistic line profile in the average 
rest- frame spectrum of type-1 and type-2 AGN. In objects with 



broad Fe lines, the contribution of the broad component is dif- 
ficult to detect if there is not enough signal to noise since its 
contribution is less than 10% above the continuum over most 
of the spectrum. 
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channel energy (keV) 

Fig. 16. Unfolded MOS and pn time averaged spectra of the two still unidentified sources (source numbers 26 and 537) where 
we detected signatures of emission line at high energies (F-test significance >95%). 
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Fig. 17. Evolution with redshift of F and Nj^J" (rest-frame) for type-1 and type-2 AGN. Error bars correspond to 90% confidence 



10. Dependence of sources spectra with 
luminosity and redshift 

We show r and N™" vs. redshift for type-1 AGN and type-2 
AGN in Fig.^l In these plots we can see the different redshift 
distributions between the type-1 and type-2 AGN in our sam- 
ple. Most detected type-2 AGN have redshifts below 1, while 
we find type-1 AGN up to a redshift of -3.5. We have applied a 
Spearman correlation test to search for evolution of F and NJ^"^ 
with redshift. We found that the correlation between F and red- 
shift is -0.22 for type-1 AGN and -0.04 for type-2 AGN. The 
significance of F being flatter at higher redshifts is 86% for 
type-1 AGN and 15% for type-2 AGN. The continuum shape of 
our sample of AGN does not seem to evolve with redshift, how- 
ever the number of AGN at high redshift (specially the number 
of type-2 AGN) is too small to give a strong conclusion. 

The same result is obtained when searching for correlation 
of N™"^ with redshift, i.e. AGN at high redshift do not seem 
to be more absorbed than local ones. The apparent scarcity of 



high redshift (z>l) low Nh sources is probably a selection ef- 
fect, since it is easier to detect highly absorbed sources at high 
redshifts. 

We have studied the dependence of F and NJ^* with the 2-10 
keV X-ray luminosity (we use the 2-10 keV luminosity because 
these values are less affected by X-ray absorption). The results 
are plotted in Fig We do not see any correlation between 
the spectral slope and the column density of our AGN with 
2-10 keV X-ray luminosity. Note in Fig. ^] that our sample 
contains 6 objects (all optically identified as type-2 AGN) that 
fall within the "standard" QS02 region, i.e., Lx> 10^"^ erg s ' 
andNjj'^ > lO^^cm^^. 

11. Unabsorbed type-2 AGN 

We have found 5 objects identified as type-2 AGN but with 
no clear evidence of X-ray absorption in their X-ray spectrum 
(~ 23%, see Table. |3}- In Fig.^]we show the unfolded spectra 
of these sources obtained with their best fit model (in all the 
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Fig. 18. Left: T vs. 2-10 keV luminosity. Right: Nj^"' vs. 2-10 keV luminosity for type-1 and type-2 AGN. In the Nj^"^ vs. 2-10 
keV plot we show the QS02 region (defined as L2-10 > 10'*"^ergs ' and Nh > 10^^ cm"^) where we have 6 QS02 candidates. 
Error bars correspond to 90% confidence 



cases a single power law). Several authors have found AGN 
with weak or no broad emission lines in their optical spectrum 
and with unabsorbed X-ray spectra (see for example Pappa et 
al., ,2001 Panessa et al., 2002. Barcons et al., 2003 , Mateos 
et al., I2OO5I Carrera et al., 120041 Corral et al. l2005l . 

One possible explanation for these results is that the signal 
to noise of the spectra of these sources is not high enough to 
detect signatures of X-ray absorption (see e.g. Mateos et al., 
I2005 >. Our current data is however of sufficient quality to de- 
tect X-ray absorption or X-ray absorption+soft excess to weak 
levels. Moreover, as we see in Fig.^]the time averaged spectra 
of these sources has enough signal to noise as to detect X-ray 
absorption with the column densities common in type-2 AGN. 
We have calculated the upper limits (at 90% confidence) to the 
X-ray absorption in these sources. They are listed in Table 
The values that we obtain are lower than the typical column 
densities found in our absorbed type-2 AGN. If these sources 
are X-ray absorbed, the values for the column density that we 
have found are consistent with arising in absorption from their 
host galaxy. 

Another possibility that might explain the lack of X-ray ab- 
sorption in these sources is X-ray spectral variability. The op- 
tical and X-ray observations of these sources have not been 
obtained simultaneously, and the X-ray absorption in these 
sources might have changed with time. Corral et al. (2005) 
studied the hypothesis of spectral variabiUty using simulta- 
neous X-ray and optical observations of the Seyfert galaxy 
Mkn993. They found the source to be X-ray unabsorbed but 
in a type 1 .9 optical. Results of a detailed study of X-ray flux 
and spectral variability on scales from months to years of the 
same sample of sources used for this work, will be presented in 
a forthcoming paper (Mateos et al. 2005b in preparation) where 
we show that X-ray variability cannot explain the lack of X-ray 
absorption in our unabsorbed type-2 AGN. 

We have checked whether these sources can be Compton- 
thick type-2 AGN. If the torus is Compton-thick to optical 
scattering, even 2-10 keV photons will not be directly seen 



and hence the direct radiation in these sources would be com- 
pletely blocked. In some cases scattered radiation (with no ap- 
parent absorption) could be the only radiation seen below 10 
keV. In Compton-thick sources, because the primary radiation 
is fainter, the equivalent width of the Kc line increases. Bassani 
et al. (1999't show a diagram of the EW versus the transmission 
parameter T, where Tis 5x/5[oui]". This diagram can be used 
to identify Compton-thick sources if the values of EW and T 
are known. We do not know the value of S\p\\\\ for our sources, 
however we have measured the EW to check in which part of 
the diagram our sources fall. To measure the EW we added a 
Gaussian line representing the iron Kcc emission to the spec- 
trum of our sources. We fixed the centroid of the line to 6.4 and 
the width to (the quality of the fits did not improve allow- 
ing the line parameters to vary). Only in source 21, there might 
be emission from iron Ka line (F-test significance ~95%). In 
the other sources there are no indications of iron Ka emission. 
The values of the EW that we have found are listed in column 
(6) of Table In most of the cases we obtained a value of the 
EW below ~1500, and therefore these sources fall outside the 
region of Compton-thick sources in the Bassani et al. J1999> 
diagram. However, to confirm our results, specially for objects 
21 and 476 with EW values above 1000 eV, we need a reliable 
measurement of the [OIII] flux. 

12. Extragalactic X-ray background 

The spectrum of the extragalactic X-ray background (XRB) 
was measured by the HEAO satellite (Marshall et al. 1980) 
from 1-50 keV. At these energies, the XRB spectrum can be re- 
produced well by an optically thin plasma of temperature ~40 
keV. At low energies, <15 keV, a good description of the data 
is obtained with a power law of r=1.4. The spectrum of the 

" Sx is the X-ray flux and 5[om] is the optical flux of the [OIII] 
/15007 emission line. [OIII] /i5007 has been frequently used as an 
isotropic indicator of the intrinsic brightness of the sources. 
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XRB is significantly flatter than the typical spectrum of AGN. 
A population of heavily absorbed AGN, predicted by synthesis 
models of the XRB (e.g. Setti & Woltjer, .1989. Comastri et al.. 



1995"; Gilli et al., |200T1 Gandhi & Fabian, l2003l Ueda et al., 
i2003i) . might account for this discrepancy. 



We have carried out a stacking of MOS/pn time averaged 
spectra of the sources that we have analysed. Details on the 



S. Mateos et al.: X-ray spectra LH 



19 



Table 7. X-ray properties of the type-2 AGN for which we did 
not found absorption in thek X-ray spectrum. 



ID 

(1) 


redshift 
(2) 


r 

(3) 


(90%prob.) 
(4) 


L2-10 
erg s"' 
(5) 


Equivalent 
width (eV) 
(6) 


6 


0.528 


1.94 ±0.20 


< 20.91 


43.07 


576 


21 


0.498 


1.77 ±0.14 


< 20.62 


43.28 


1462 


39 


0.711 


1.79 ±0.08 


< 20.66 


43.44 


292 


427 


0.696 


1.40 ±0.20 


< 21.08 


42.79 


249 


476 


0.607 


2.07 ±0.17 


< 20.22 


42.66 


1393 



Columns are as follows: (1) Source X-ray identification number; (2) 
redshift; (3) T from best fit model (for all sources the best fit model 
was a single power law); (4) upper limit in intrinsic (rest-frame) X-ray 
absorption (90% confidence); (5) logarithm of the 2-10 keV luminos- 
ity; (6) rest-frame equivalent width of an emission line centred at 6.4 
keV with <t=0 (the value of the EW was obtained using a Gaussian to 
fit the emission line) 

analysis are given in AppendixIXl The goal of this analysis was 
to compare the integrated emission of our sources, with an av- 
erage spectral shape of ~1.92 (see Sec. l9.1> . with the spectrum 
of the XRB in the 2-7 keV energy band. 

We fitted MOS and pn stacked spectra with xspec using 
a simple power law model. Then we divided the measured 
XRB intensities (in units of keV^keV"' cm"^ s"' at 1 keV) 
by the total solid angle covered by XMM-Newton (0.4 deg^ or 
1.2185x10 sr) to obtain the resolved fraction of XRB in the 
area surveyed. 

Fig.|20|shows the total extragalactic XRB spectrum as mea- 
sured by the HEAO-1 mission (solid line) but renormalised to 
the 2-8 keV intensity observed by De Luca & Molendi (2004i. 
The points show the 2-10 keV XRB spectrum seen by MOS 
(diamonds) and pn (stars). 

We found that the 2-7 keV XRB resolved by our 
sources was best fitted with a power law of F = 
1.59 ± 0.03 and N=4.78 + 0.15keV2keV-i cm^^ g-' sr"' (at 
1 keV) for MOS data and Y = 1.54 + 0.04 and 
N=4.38 + 0.14keV2keV-icm-2s-' sr-i (at 1 keV) for pn 
data. The integrated contribution of our sources is indeed 
harder than the spectrum of the brightest AGN, but still 
softer than the XRB at these energies. Since our sources 
have been selected in the 0.2-12 keV band, we probably 
missed faint absorbed sources. Worsley et al. ( i2004i l did in- 
clude fainter sources than the present study, and hence they 
reached a significantly higher integrated emission at 1 keV 
(N = 1 1 + 0.5 keV^ keV"' cm"^ s"' sr"'). A further component 
not detectable by XMM-Newton arising above 5 keV might also 
be present, as argued by Worsley et al. J2004t . 

13. Discussion and Conclusions 

We have caiTied out a detailed study of the X-ray spectra 
of a sample of 123 objects detected with XMM-Newton in a 
deep observation in the Lockman Hole field. The EPIC spec- 
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Fig. 20. Spectrum of the extragalactic X-ray background as 
measured by the HEAO satellite but renormalised to the 2-8 
keV intensity observed by De Luca & Mondeli (2004) (solid 
line). The points show the 2-7 keV stacked spectra of the 
sources that we have analysed (M14-M2 (diamonds) and pn 
(stars)). 

tra of these sources have all more than 500 background sub- 
tracted counts (MOS-npn) in the energy interval form 0.2-12 
keV. Thanks to the good signal to noise of the data we could 
study in detail the different spectral components that contribute 
to the X-ray emission of AGN in the above energy interval. 

13.1. X-ray continuum shape and intrinsic absorption 

The 0.2-12 keV spectra of many AGN cannot be well repro- 
duced with a single power law model. We found that this model 
is the best fit model for only ~53% of the sample. The average 
continuum shape of our sources appears to harden at fainter 
0.5-2 keV fluxes (see e.g. Giacconi et al., 2001,. Mainieri et 
al.. l2002l Mateos et al.,|^5), but there are some indications 
that the same effect is seen for objects detected in the 2-10 keV 
band (Tozzi et al., .200L Streblyanska et aL. l2004b . Two hy- 
pothesis have been suggested to explain it: X-ray absorption is 
more important at fainter fluxes or there exists a population of 
faint sources with intrinsically flatter continuum. 

We confirm that this apparent effect is due to absorption. 
Both because it is not seen in the 2-10 keV band, and because 
it disappears when absorption is taken into account. We do not 
see any evidence for a population of intrinsically harder sources 
at faint fluxes. However, if sources with larger column densi- 
ties were found, then we should see the hardening of (F) in 
the 2-10 keV band. Using a single power law model. X-ray 
absorbed sources had the faintest 0.5-2 keV fluxes, but they 
had the same distribution of 2-10 keV fluxes than unabsorbed 
sources. Because of this, the average F obtained using a single 
power law model as a function of the 2-10 keV flux was mea- 
sured to be harder (F ~ 1.7) than the typical value of ~ 1.9 
found in unabsorbed AGN. 

In the hard band we found no dependence of X-ray absorp- 
tion or fraction of absorbed sources with flux. If observed (not 
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corrected for absorptioti) 0.5-2 keV fluxes are used, theti the 
column density and the fraction of absorbed sources increase 
significantly at fainter fluxes. However, using unabsorbed 0.5- 
2 keV fluxes (i.e. fluxes corrected for the eff'ect of absorption) 
the correlations are not observed and we obtain the same results 
than in the 2-10 keV band. Therefore to study the dependence 
of r with the X-ray flux it is important to first correct the fluxes 
for the effect of X-ray absorption. Another interesting result 
from our analysis is that allowing the objects to be X-ray ab- 
sorbed we found a value of (F) ~ 1.9 at all 0.5-2 keV and 2-10 
keV fluxes but with significant scatter in the points. 

We have carried out a detailed study of the X-ray spectra of 
each individual object. We have searched for soft excess emis- 
sion and for signatures of reflection at high energies (Compton 
reflection). We also studied the presence of Fe Ka emission. 
Using for each source the value of F obtained after including all 
these components to the fitting model we found that the scat- 
ter in (F) is much smaller, and therefore is mostly due to the 
presence of other spectral components in the X-ray emission. 

However we still expect some intrinsic scatter in F. To cal- 
culate the intrinsic dispersion in the continuum shape of our 
objects, we assumed that the distribution of F could be well 
represented with a Gaussian. Under this hypothesis, we found 
our sources to have an average spectral slope of ~ 1.92 with an 
intrinsic dispersion of ~0.28. 

13.2. Soft excess emission 

Soft excess was detected in the time averaged spectra of 18 
sources. However only in 9 objects (4 type-1 AGN, 4 type- 
2 AGN and 1 unidentified source) we could fit the spectral 
signatures with a black body model (a Raymond Smith gave 
an equally good fit). We found the average temperature of the 
black body to be 0.09+0.01 keV for type-1 AGN andO.26+0.08 
keV for type-2 AGN. The average 0.5-2 keV luminosities of the 
black body were found to be (in log units) 43.42 ± 0.43 erg s"' 
in type-1 AGN and 44.11 + 0.44 erg s ' in type-2 AGN. The 
0.5-2 keV luminosities of the soft excess component do not 
differ significantly between type-1 and type-2 AGN, but our 
results seem to indicate that the black body temperatures are 
slightly higher in type-2 AGN than in type-1 AGN. This might 
be due to a higher contribution from scattering in type-2 AGN. 
However due to the small number of AGN with soft excess 
analysed we can not reach any conclusion. The temperatures 
of the black body are in most cases >60 eV, and therefore the 
origin of the soft excess component in these sources cannot 
be explained as thermal emission from the accretion disc only. 
Componization of cool disc photons by hot electrons surround- 
ing the accretion disc might be an alternative explanation. 

In 9 sources (including 1 type-1 AGN, 3 type-2 AGN and 
5 sources unidentified) the black body model could not fit 
the spectral signatures of the soft excess emission. For these 
sources we obtained a good fit with a scattering or partial cov- 
ering model (two power laws with the same spectral index but 
different absorptions). The average covering fraction of the ab- 
sorber was found to be 0.82 ± 0.06, which means that the scat- 
tering fraction in these sources is rather large (18 ± 6%). 



13.3. Reprocessed components 

We found one object in our sample with a flattening in the spec- 
tral slope at high energies. The source is still unidentified so we 
used a model with two power laws to fit its X-ray emission. 

Signatures of an emission line at high energies were found 
(F-test significance >95%) in 8 sources (1 type-1 AGN, 5 type- 
2 AGN and 2 unidentified sources). Although in some sources 
the profiles of the lines did not appear to be symmetrical, the 
signal to noise of the spectra was not high enough to use more 
physical models. Therefore we fitted the lines in all the sources 
with a Gaussian model. For most AGN we found the centroids 
of the lines to be consistent (within the error bars) with being 
Fe Ka at an energy of ~6.4 keV (rest frame). Only in one ob- 
ject we found the line centroid to be sUghtly higher than the 
expected value for neutral iron. In this source the line might 
arise as reflection in an ionised accretion disc. 

In two objects we found significant widths for the emission 
lines. In another two sources the centroids of the lines were 
found at energies lower than 6.4 keV. In all these sources the 
emission lines might have relativistic fine profiles, and there- 
fore they might have been emitted in the inner regions of the 
black hole accretion disc. 

13.4. X-ray spectra of AGN 

We found the best fit average spectral slope to be 1.89 ± 0.03 
in type-1 AGN and 1.71 + 0.03 in type-2 AGN. These values 
seem to indicate that type-2 AGN have harder spectral slopes 
than type-1 AGN. However assuming that the values of F that 
we have obtained for each individual object follow a Gaussian 
distribution, and allowing the spectral slope to have intrinsic 
dispersion, the significance of type-2 AGN being harder than 
type- 1 AGN is of 1 .62cr. This small difference in (F) for type- 1 
and type-2 AGN might be due to the fact that if the signatures 
of absorption are not very significant, the detected values of 
Nh will tend to be lower than the real ones and then, the fitted 
F will be flatter. This effect will be more important for type-2 
AGN where we expect more sources to be absorbed. 

X-ray absorbed objects were found among type-1 AGN 
(~ 10%) and type-2 AGN (~ 77%). We found the fraction of 
absorbed objects in type- 1 and type-2 AGN to be different with 
a significance of > 99.99%. The distribution of absorbing col- 
umn densities also suggest that type-1 AGN are less absorbed 
than type-2 AGN. A comparison with a KS test gave a signifi- 
cance of these distributions being different at 92% confidence. 

We did not see a dependence of the AGN continuum shape 
with the X-ray luminosity or redshift. We found the same re- 
sults for the column density. 

13.5. Unabsorbed type-2 AGN 

We studied in more detail the X-ray emission of the 5 type-2 
AGN with unabsorbed MOS and pn X-ray spectra. AH these 
sources have spectra with enough signal to noise, hence we 
should have been able to detect signatures of X-ray absorption 
if they are present. We have calculated the upper limit (90%) in 
the column density. We found that in aU cases the values for the 
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column density are significantly lower than the typical values 
found in the absorbed type-2 AGN. If there is X-ray absorption 
in these sources, the low values of the column densities that we 
find could be explained as arising from the host galaxy. We ar- 
gue that spectral variabiUty is unlikely to be at the heart of the 
apparent X-ray/optical mismatch. 

Finally we do not find compelling evidence that these 
sources are Compton-thick, although the [OIII] flux needs to 
be measured to reach a firm conclusion. 
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Columns are as follows: (1) XMM-Newton identification number; (2) /?05AT identification number; (3) Right ascension (J2000); (4) Declination (J2000); (5) Optical class 
from optical spectroscopy; (6) Source redshift; (7) Best fit model of the X-ray spectrum of each object: SPL: single power law, APL: absorbed power law, SE: Soft-excess 
emission, CAPL: absorption with partial covering absorber, 2SPL: two power laws; (8) Slope of the broad band continuum emission; (9) Logarithm of the hydrogen colunm 
density of the X-ray absorber (observed Nh if the source is stUl unidentified or rest-frame absorption if the source is identified; (10) and (1 1) 0.5-2 and 2-10 keV flux in units 
of 10"'"* erg cm^ s ' obtained from the best fit model; (12) and (13) logarithm of the 0.5-2 and 2-10 keV luminosities (corrected for absorption) for the identified sources 
obtained from the best fit model. 

The source coordinates RA and DEC, are not the X-ray positions from the Lockman Hole catalogue (H. Brunner et al. 2005, in preparation), but the centres of the 
regions used to extract the spectra of the objects. 
* Source photometric redshift. 
** BAL QSO 
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Appendix A: Stacking of spectra of Lockman Hole 
sources 

To calculate the contribution to the Cosmic X-ray background 
from our sources, MOS and pn spectra have been stacked. 

We have kept MOS and pn data separately because of the 
different instrumental responses and to provide support to the 
results. In addition, as it was explained in Sec.|6l MOS and pn 
time averaged spectra were not necessary built using the same 
set of observations, hence the normalisations of the XRB spec- 
trum obtained with MOS and pn might be sUghtly different. 

At the time of this study there were calibration uncertainties 
between MOS and pn data at energies below ~1 keV (pn gives 
higher fluxes below 0.7 keV by 10-15% with respect to MOS, 
see Kirsh et al.'2004'). We also found our data to be rather un- 
certain at energies > 7keV. Because of the rapid decrease in the 
effective area of the X-ray detectors at energies > 5keV the sig- 
nal to noise of the data becomes very low. In addition we know 
that the particle background is very important at high energies. 
Therefore we restricted our analysis to the 2-7 keV energy in- 
terval where we know that our results will not be affected by 
calibration problems or inaccurate background subtraction. 

MOS and pn stacked spectra were obtained with the fol- 
lowing procedure: suppose we want to add two spectra hav- 
ing exposure times tl and t2, backscales bl and b2 and re- 
sponse matrices rspl and rsp2, and that the corresponding 
background extraction regions have exposure times tbl and 
tb2 and backscales bbl and bb2: 

1. Add spectra: The total backscale will be bl+b2'- and the 
total exposure time Ctlxbl+t2xb2)/(bl+b2). 

Note that we want to extract the XRB from the regions where 
we have the sources, and therefore, we are accumulating counts in an 
increasing solid angle. 
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2. Add background files: The total backscale 
will be bbl+bb2 and the total exposure time 
(tblxbbl+tb2xbb2)/(bbl+bb2). 

3. Add response matrices. The combined response matrix will 
be Ctlxbl)xrspl+Ct2xb2)xrsp2/Ctlxbl+t2xb2). 

This procedure can be trivially extended to stack any number 
of sources. 



